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The local Zjv quantized Berry phase for the antiferromagnetic Heisenberg spin model is 

formulated. This quantity, which is a generalization of the local Z 2 Berry phase for SU(2) symmetry, 
has a direct correspondence to the number of singlet pairs spanning on a particular bond, and is 
effective as a tool to characterize and classify the various symmetry protected topological phases of 
one-dimensional S\J{N) spin systems. We extend the path-integral quantum Monte Carlo method 
for the Z2 Berry phase in order to calculate the Zjv Berry phase numerically. We demonstrate 
our method by calculating the Z4 Berry phase for the bond-alternating SU(4) antiferromagnetic 
Heisenberg chain, which is represented by the Young diagram of four columns. 

PACS numbers: 02.70.Ss, 03.65.Vf, 05.30.Rt, 75.10.Pq 


The ground state of S' = 1 antiferromagnetic Heisen¬ 
berg (AFH) chain is unique and has a finite energy gap, 
while that of the classical AFH chain is continuously de¬ 
generated and has no energy gap. The former phase, 
called Haldane phase [T], is understood by the valence 
bond solid (VBS) picture [2] and characterized by the 
hidden string order [3]. The Haldane phase emerges in 
other one-dimensional systems, such as the spin ladder 
systems HIS], and has been studied for a long time by 
means of several topological order parameters, such as 
the string order parameter HISllTl, the twist order pa¬ 
rameter [8], and the local Z 2 Berry phase [9]. 

Recently the Haldane phase has been reexamined from 
a more general viewpoint, called symmetry protected 
topological (SPT) phase [TU]. In one-dimensional sys¬ 
tems with short range interaction, all topological phases 
are adiabatically connected to a trivial state (product 
state) by changing the parameters if one does not im¬ 
pose any symmetries to the Hamiltonian; some symme¬ 
tries are required to make topological phases nontrivial 
and distinguish them from the trivial state. For exam¬ 
ple, the Haldane phase of the S' = 1 AFH chain is pro¬ 
tected by any one of the following three symmetries m- 
time reversal symmetry, bond-centered inversion symme¬ 
try, or the dihedral group of tt rotation around the Sx, 
Sy, and Sz spin axes. As a promising tool for characteriz¬ 
ing such topological phases, the entanglement spectrum 
of the ground state has been proposed HH. This quan¬ 
tity is the spectrum of the reduced density matrix, which 
is obtained by tracing out the degrees of freedom in one 
part of the whole system from the original density ma¬ 
trix. In the Haldane state of the one-dimensional AFH 
S' = 1 chain, the lowest eigenvalue in the entanglement 
spectrum is doubly degenerated, while that of the trivial 
state is unique m- 

In the meanwhile, as a candidate that exhibits novel 
SPT phases, spin models with generalized S\J{N) sym¬ 
metries have attracted attentions in years. In general, the 


SPT phases of one-dimensional spin systems with some 
symmetry group G are classified by the second cohomol¬ 
ogy group U) of G [T5HT5] . It is pointed out that 

a one-dimensional SU(A^) spin system can have A^—1 dis¬ 
tinct SPT phases besides the trivial phase without any 
additional symmetries m. though the SU(2) S' = 1 spin 
chain has only one nontrivial topological state, the Hal¬ 
dane state. If we cut a bond in the SU(A^) spin chain, 
SU(A^) edge states emerge at both ends. The SPT phases 
of the S\J{N) spin chain can be classified by the number 
of boxes in the Young diagram representing these edge 
states. In other words, if we have a means to count the 
number of boxes directly, the SPT phases can be identi¬ 
fied explicitly. For example, the degeneracy of the lowest 
eigenvalues in the entanglement spectrum is the same as 
the dimension of the representation of the edge state, 
which helps us to estimate this number. 

So far, a number of SU(A^) models have been pro¬ 
posed to study the SPT phases. Morimoto and his 
co-workers [TB] constructed a matrix product state rep¬ 
resenting an S\J{N) VBS state from the viewpoint of 
SPT and the corresponding Afheck-Kennedy-Lieb-Tasaki 
(AKLT) Hamiltonian. They also calculated the string or¬ 
der parameter [17] of the SU(3) AKLT model by the den¬ 
sity matrix renormalization group method. Geraedts and 
Motrunich [TB] proposed another exactly solvable SU(V) 
model that has N distinct SPT phases as the ground 
state by generalizing the cluster Ising model [HI. Quest 
for S\J{N) SPT states in the real world is also in progress. 
Ultracold alkaline-earth atom systems, such as ®^Sr, 
^^^Yb, or ^^^Yb atoms, in an optical lattice are described 
by a two-orbital SU(V) Hubbard model or an SU(V) 
AFH model in the Mott insulating phase [20] • Nonne 
and his co-workers studied the former model analytically 
and numerically, and found SU(V) SPT phase [2T] . 

In this Letter, we introduce the local Berry 
phase [22] to the SU(V) AFH model and show that this 
quantity can detect the number of SU(V) singlet pairs 
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on a particular bond. This number corresponds directly 
to the number of boxes in the Young diagram of the 
edge states, and thus enables us to distinguish several 
SPT phases. We then extend the path-integral quantum 
Monte Carlo (PIQMC) method for calculating the local 
Z 2 Berry phase [53] to the present Z^v case. To demon¬ 
strate the power of the Berry phase and our PIQMC 
method, we perform a simulation of the bond-alternating 
SU(4) AFH chain, which is represented by a Young di¬ 
agram of four columns, and confirm that the local Z 4 
Berry phase indeed identifies the nontrivial SPT phases 
correctly. 

The SU(fV) antiferromagnetic spin model can be gen¬ 
erally written as 


the twist vector 0 , whose elements are defined as 

l<a<N 

0 a = N. ^ ’ 

By using (j), we “twist” the dimer Hamiltonian as 

= - — -(/3^|. (4) 

a,/3 

The corresponding ground state is also twisted as 

Next, we param¬ 
eterize {0a} by introducing a “time” parameter t as fol¬ 
lows: 


N 

{i,j) a,/3=l 


( 1 ) 

and 


9i — t 


( 5 ) 


where (i, j) represents a nearest-neighbor bond connect¬ 
ing sites i and j, and the outer summation runs over all 
the nearest-neighbor bonds. Sp denotes a generator of 
SU(A^) algebra satisfying the following commutation re¬ 
lation: [S^{i),Sl{{j)] = In this 

Letter, we consider bipartite lattices and adopt the rep¬ 
resentation of the algebra depicted as a Young diagram 
of one row and M columns for one sublattice A and the 
one depicted as one of (N — I) rows and M columns for 
the other sublattice B |53H53] . For M = I, the model 
can be written in terms of an Wcolor particle basis as 

^ = - E E E {m I > (2) 

i^A j^B cxi^otj 

where \ai) = ^ |vacuum) and ^ is the fermion 

creation operator of color a at site i, and |dj) = 
Ca j 11/3 4.1 I’^acuum) = (-1)“ 0 / 3/0 |/3j) ■ Physically, |a) 
and 15) mean an a particle and hole, respectively |24j . 
When Jij > 0, the ground state of the dimer Hamilto¬ 
nian has a unique ground state |^) = 
with energy E = —\ and — I)-fold degenerated ex¬ 
cited states with E = 0. We call this ground state as the 
SU(A^) singlet. 

The ai = 02 subspace of the Hamiltonian of SU(A^) 
AFH dimer, which is represented by an x matrix 
with all the elements being — I/A and still includes the 
SU(A) singlet state, is the same as the Hamiltonian of a 
hopping particle on A fully connected sites. Since the 
ground state of the latter is characterized by the Zjv 
Berry phase |55|, so is the SU(A) singlet bond. In the 
following, we will introduce the Z/v Berry phase 7 for the 
present SU(A) AFH model as in Ref. |22| and show ex¬ 
plicitly that the value of the Berry phase of the SU(A) 
singlet is given by 7 = —27r/A mod 27r. First, we intro¬ 
duce (A — 1) parameters 9a with a = I,2,...,A—I and 


t 0 < t < 2 tt/N 
2'k/N < t < 2 tt 


( 6 ) 


for 1 < a < A. By using this parameterization, we define 
the Berry curvature Aa{(l>) as 


Aa{$) 




( 7 ) 


and the Berry phase 7 as 


7 = d$- A{$) = i J dt'^^^Aa{(j){t)). 


( 8 ) 


It should be noted that (j>a{2Tr) = 0 mod 27r and 
'H{(j){2'K)) — make the parameter path a 

closed loop and accordingly the Berry phase becomes 
gauge invariant. Finally, the Berry curvature and the 
Berry phase of the SU(A) dimer can be evaluated ex¬ 
plicitly as 


Aa{$) = 





and 


i 

A 


( 9 ) 



= 27r 


A- I 
A 


( 10 ) 


respectively, where we used relations 4>n = 0 and 
0/3 (27r) - 0/3(0) = 27rJ/3,i. 

For M > I, a spin can be regarded as a symmet¬ 
ric summation of M one-column SU(A) subspins as an 
SU(2) spin with S' > 1/2 can be decomposed into 2S 
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subspins with S = 1/2. Thus, the ground state of M- 
column SU(fV) dimer is the symmetric superposition of 
M\ product states of M singlet pairs as 

, M N (11) 

where / is the index of subspins and P denotes a per¬ 
mutation of M distinct elements. A short explicit cal¬ 
culation gives that the cross term in the Berry con¬ 
nection, , becomes —iM/N ('Ppj'I'p/) [57]. 

The Berry connection and the Berry phase are obtained 
in the same way as the M = 1 case as Aa{4)) = 
(^'((/)|9a I '!'(</)) = —iM/N and 7 = 2ttM{N — 1)/A^, 
respectively. 

The value of the local Berry phase on bond b thus 
gives the number of singlet pairs on the bond as n{b) = 
'){h)N/{2 tt{N — 1)). If we cut bond b, n{b) singlets are 
removed and 2n{b) subspins are left as the edge state. 
The edge state on the site belonging to sublattice A is 
depicted as the Young diagram of 1 x n(b) boxes. For 
one-dimensional spin systems, the number of boxes char¬ 
acterizes SPT phases [15] , and so does the local Zjv Berry 
phase. 

In the following, we will describe how to calculate 
the Berry connection and the Berry phase by means of 
the PIQMC method. The Berry connection, A(t) = 
('i/(t)jdtj'i/’(t)}, is the coefficient of the first order term 
of the series expansion of the inner product between two 
normalized ground states: 

+ St)} = 1-1- StA{t) + 0{5t^). (12) 


The ground state \'4’{t)) of the Hamiltonian 'H(t) can 
be obtained by the projection method as |^'(t)) = 
lim^_).oo Y(t, ,5) exp (—/3'H(t)/2) \(j )), where (3 is the pro¬ 
jection parameter, N(t, /3) G K a normalization fac¬ 
tor, and \(j)) some initial state nonorthogonal to 
the ground state. We will omit the “lim” sym¬ 
bol for simplicity from now. The inner prod¬ 
uct of two ground states ( 12 ) is also represented 
by the projection as + St)) = N{t)N{t + 

St) {4>\ exp {—(3V.{t)/2) exp + St)/2) \(j)). 

Next, as one expands Tr exp(—/3H) in the ordinary 
PIQMC, we expand {cj)\ exp(... )|^) by path integrals as 


{ip{t)\tp{t + St)) = N{t)N{t+St) ^ W{c] t,t+St), (13) 


where W{c\t,t + St) is the weight of a worldline config¬ 
uration, c. By using the identities, N{t)"^ X)c ^('-j () = 

m)\m) = 1 and N{t)A-^N{t) = ^lV(t)V2, Eq. 
can be expanded with respect to St up to the first order 


r 



FIG. 1. (Color online) Examples of worldline configuration 
(left panel) and loop configuration (right panel) of the SU(A) 
AFH model. Color of worldlines denotes spin state; black line 
(A), red dashed line (1), and blue dotted line (2). 


as 


{'ii){t)\ii){t + St)) 

'5^Ec(^tL - 5tu)lF(c;tu,tL) 


~ 1 -b 


= 1 -I- 


2 EcW{c;t,t) 

(dtj - (dtj 


tL = tu=t (14) 


where (O) = 0(c)VF(c, t, t)/lF(c, t, t). Combin¬ 

ing Eq. (14) with Eq. (12) we obtain a generic form for 
the Berry connection. 


A{t) = -{{d,,)-{dt,)). 


(15) 


Note that since we do not use any knowledge of any spe¬ 
cific models so far, this formula is general and valid for 
any models. 

To advance the procedure further, we now need an ex¬ 
plicit representation for (d) based on the weight function 
W(c) of the model under consideration. The expecta¬ 
tion of derivative with respective to t can be evaluated in 
the same way as the energy evaluator, ("H) = — ( 9 / 3 ) = 
— (to) //3, where to is the number of operators |55|. For 
the present SU(A) AFH model (|^, the evaluator of the 
Berry connection is eventually written as 


N 


Me; <?) = -X (16) 


3=1 


where N^^'^ is the number of the operators, lad) {I3j3\, on 
the twisted bond at the imaginary time r e [0, /3/2] (t e 
[P/2,13]). The left panel of Fig. shows an example of 
worldline configuration of four site system. Black lines, 
red dashed lines, and blue dotted lines denote spin state 
N,l, and 2, respectively. When we twist bond (1,2), 
= N// = 1 and the other N/^^s and N^/^s are 

zero in this configuration. Note that the twist pertur¬ 
bation makes the weights complex-valued, so we should 
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use the reweighting method to perform Monte Carlo sim¬ 
ulation as (A) = {AS)q / {S)q, where (•)g is the ex¬ 
pectation value under the absolute weight system with 
Wq = |VC|, and S = Wl\yV\ is the phase function written 
as S'(c;0) = ria,^ 6 xp[—*(</)„ — (j)i3)Nap] for the present 
model. Berry connections for all parameters can be cal¬ 
culated from one Monte Carlo simulation, and the Berry 
phase can be obtained by integrating these numerically. 

For updating the worldline configuration, we employ 
the loop cluster algorithm [281130] . In the loop cluster 
algorithm, one can represent the Berry connection as a 
function of loop configuration (improved estimator) in¬ 
stead of the worldline configuration [23] • It should be 
noted that in the present projector Monte Carlo, the 
states at T = 0 and /3 are fixed to some reference state. 


which break edge loops and prevent these from flip¬ 
ping. We choose \(j)) as the “Neel” state, as 

this choice makes the improved estimators simple be¬ 
cause oi (j}N — 0 for any t. In the improved estimator, 
the Berry connection and the phase of the weight are de¬ 
composed into the contributions from each closed loop 
as 2iAa = and S = expi0Q(^)?^-i-(^), 

respectively, where a{£) is the state (color) of loop £ and 
n± = (IVl.t - ^L.;) ± (IVu,t - ^u.;), where these N’s 
are the numbers of the jumps above (f) and below (),) 
the operator at r < /3/2 (L) and t > 13/2 (U) on the 
twisted bond. By tracing out the spin state, we obtain 
the expression of the improved estimator for the Berry 
phase as 


7 =* rdt^ ■!{${£))= 

Jo 


dt 


dt 


? EE 




^—in+{£')(p. 


dt 


N 



where the loop indices £ and £' run over all closed loops. 
For N = 2, this estimator is reduced to the one derived 
in the past work [23] . A typical loop configuration for the 
SU(A^) AFH model is shown in the right panel of Fig. 

If we twist bond (1,2), n± of the thick loop are both 1, 
and those of the other loops are all 0 in this case. 

We can extend the above procedure straightforwardly 
to higher-order terms in Eq. (12|, which yields estima¬ 
tors for other useful quantities. From the second or¬ 
der, for example, the expressions for the susceptibil¬ 
ity of fidelity [3I| xf = + St)f /{6t)‘^ 

and the Berry curvature flap = daAp — dpAa are 
derived as xf = {d^di/} - (5u) (^l) and flap = 

{{dapjdpx) {da,v) {dpx)) - (a o /3), respectively. We 
note that more specialized QMC estimators have been 
proposed for the susceptibility of the fidelity | 33 | and for 
the Berry curvature [33| , which are more accurate for the 
specific models. 

In order to demonstrate the ability of the local Zjv 
Berry phase, we applied our PIQMC method to the 4- 
column SU(4) AFH model {N = 4 and M = 4) on a 
bond alternating chain [24] . The Hamiltonian is written 
as 


L 

H(5)=5](1-(-I)*5)Hm+i, (18) 

i=l 

where is the two-body Hamiltonian on a bond con¬ 

necting j-th and {i -I- l)-th sites, L the length of chain, 
and the periodic boundary condition, S^{i + L) = Sp{i), 
is imposed. In the 5 = 1 limit, the ground state is obvi¬ 
ously the product state of decoupled L/2 dimers so that 
the number of singlet pairs on a 1 -I- 5 bond, n{+), is 4 


and one on a 1 — J bond, n(—), is 0. In the 5 = —1 limit, 
on the other hand, n{+) and n(—) are 0 and 4, respec¬ 
tively. As 5 is increased from —1 to 1, it is expected that 
quantum phase transitions between different topological 
phases occur successively. We calculated the Berry phase 
7 by numerically integration of the Berry connections 
A{t) evaluated by the present PIQMC method at twist 
parameters t = 27rfc/32 with k = 0,1,..., 31. For the 
chain lattice, n+(I') in Eq. ( [l7| ) is almost always zero be¬ 
cause n^(£) is nothing but the winding number of the 
loop. As a result, the phase factor of the weight, S, stays 
almost always 1, and allows us to perform a precise sim¬ 
ulation without introducing any additional technique for 
solving the sign problem, such as the meron cluster algo¬ 
rithm 133] • We also calculated the “staggered” sus¬ 
ceptibility, X) the response function against the field that 
favors the Neel state, i.e., (^i^A l«)i \d)j, by using 
the standard loop algorithm. The projection parameter 
for calculating the Berry phase and the inverse tempera¬ 
ture for the susceptibility are both chosen as (3 = 2L. 

Figure [^ shows the simulation results for the L = 64 
chain. It is clear that the staggered susceptibility tends 
to diverge at <5 ~ —0.42, —0.14, 0.14, and 0.42, suggest¬ 
ing quantum phase transitions between different phases. 
Although such successive phase transitions have already 
been observed in the SU(2) spin system with M = 4, 
i.e., S = 2 bond-alternating Heisenberg chain |5], 5- 
dependence of the Berry phase in Fig. [^elucidates that 
the nature of the intermediated phases are completely 
different from that in the SU(2) case. On the 1 — <5 bond 
(red symbols in Fig. [^, the Berry phase shows a step¬ 
like behavior and the height of each step is —3tt/2 (or 
tt/2 mod 27r), which corresponds to the decrease of the 
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FIG. 2. (Color online) 5-dependence of the Berry phase (left 
axis) and the staggered susceptibility (right axis) of the bond- 
alternating 4-column SU(4) AFH chain with L = 64. Red 
symbols forming climbing stairs (blue forming down stairs) 
denote the Berry phase 7 on a 1 — 5 (1 -f 5) bond. Black 
symbols denote the staggered susceptibility x- Both the pro¬ 
jection parameter and the inverse temperature are /3 = 2L. 
The error bar is smaller than the symbol size. 


number of singlet pairs on the bond, while on a l-|-(5 bond 
(blue symbols), the Berry phase decreases by 7r/2 (or in¬ 
creases by 37r/2) at each critical point. Since the Berry 
phase is quantized by symmetry as long as the energy 
gap is opened, two gapped states with different Berry 
phase are topologically different [3]. Thus, we conclude 
that the present model exhibits four different topological 
phases with 7 = 0, 7r/2 , tt, and 37r/2, respectively. This 
is in sharp contract with the S = 2 SU(2) case, where 
only one nontrivial topological phase exists besides the 
trivial one. 

In summary, we defined the local Zjv Berry phase for 
the SU(A^) antiferromagnetic Heisenberg model, which 
can be directly related to the number of the singlet pairs 
on a bond. We also extended the path-integral quan¬ 
tum Monte Carlo method for local Z 2 Berry phase of 
SU(2) Heisenberg model to general models, and derived 
the evaluator of the Zjy Berry phase for the SU(A^) anti¬ 
ferromagnetic Heisenberg model. For demonstration, we 
calculated the Berry phase of the SU(4) antiferromag¬ 
netic Heisenberg model for the representation depicted 
as Young diagrams of 4 columns on bond alternating 
chains. Note that in the M = 4 SU(4) case the de¬ 
grees of freedom on one spin is 35, which corresponds 
to 5.13 S = 1/2 spins in the SU(2) case. The rapid in¬ 
crease of the size of the Hilbert space makes it impossible 
to perform simulations based on wave functions, such as 
the exact diagonalization and the DMRG method. Al¬ 
though the classical quantities such as the susceptibility 
can not identify different topological phases, the Berry 
phase successfully reveals that the present model indeed 


has N distinct SPT phases. 

The simulation code used in the present study has been 
developed based on the ALPS library [351ISS]- The au¬ 
thors thank the Supercomputer Center, the Institute for 
Solid State Physics, the University of Tokyo for the fa¬ 
cilities and the use of the SGI Altix ICE 8400EX. The 
authors acknowledge support by the Computational Ma¬ 
terials Science initiative (CMSI) in the HPCI Strategic 
Programs for Innovative Research (SPIRE) from MEXT, 
Japan. S.T. acknowledges the support by KAKENHI 
(No. 23540438, 26400384) from JSPS. 


[1] F. D. M. Haldane, Phys. Lett. A 93, 464 (1983). 

[2] I. Affleck, T. Kennedy, E. H. Lieb, and H. Tasaki, Phys. 
Rev. Lett. 59, 799 (1987). 

[3] M. den Nijs and K. Rommelse, Phys. Rev. B 40, 4709 
(1989). 

[4] H.-J. Mikeska and A. K. Kolezhuk, in Quantum Mag¬ 
netism, edited by U. Schollwock, J. Richter, D. J. Par¬ 
nell, and R. F. Bishop (Springer-Verlag, Berlin, 2004) 
pp. 1-83. 

151 T. Hakobyan, J. H. Hetherington, and M. Roger, Phys. 
Rev. B 63, 144433 (2001). 

[6] M. Oshikawa, J. Phys. Condens. Matter 4, 7469 (1992). 

[7] S. Todo, M. Matsumoto, C. Yasuda, and H. Takayama, 
Phys. Rev. B 64, 224412 (2001). 

[8] M. Nakamura and S. Todo, Phys. Rev. Lett. 89, 077204 

( 2002 ). 

[9] Y. Hatsugai, J. Phys. Soc. Jpn. 75, 123601 (2006). 

[10] Z.-C. Gu and X.-G. Wen, Phys. Rev. B 80, 155131 
(2009). 

[11] F. Pollmann, A. M. Turner, E. Berg, and M. Oshikawa, 
Phys. Rev. B 81, 064439 (2010). 

[12] H. Li and F. D. M. Haldane, Phys. Rev. Lett. 101, 010504 
(2008). 

]13] X. Ghen, Z.-G. Gu, and X.-G. Wen, Phys. Rev. B 83, 
035107 (2011). 

[14] X. Ghen, Z.-G. Gu, and X.-G. Wen, Phys. Rev. B 84, 
235128 (2011). 

[15] K. Duivenvoorden and T. Quella, Phys. Rev. B 87, 
125145 (2013). 

[16] T. Morimoto, H. Ueda, T. Momoi, and A. Furusaki, 
Phys. Rev. B 90, 235111 (2014). 

[17] K. Duivenvoorden and T. Quella, Phys. Rev. B 86, 
235142 (2012). 

[18] S. D. Geraedts and O. 1. Motrunich, arXiv:1410.1580vl 

[19] J. K. Pachos and M. B. Plenio, Phys. Rev. Lett. 93, 
056402 (2004). 

[20] A. V. Gorshkov, M. Hermele, V. Gurarie, G. Xu, P. S. 
Julienne, J. Ye, P. Zoller, E. Dernier, M. D. Lukin, and 
A. M. Rey, Nat. Phys. 6, 289 (2010). 

[21] H. Nonne, M. Moliner, S. Capponi, P. Lecheminant, and 
K. Totsuka, Europhys. Lett. 102, 37008 (2013). 

[22] Y. Hatsugai and 1. Maruyama, Europhys. Lett. 95, 20003 

( 2011 ). 

[23] Y. Motoyama and S. Todo, Phys. Rev. E 87, 021301(R) 
(2013). 

[24] I. Affleck, Phys. Rev. Lett. 54, 966 (1985). 

[25] N. Read and S. Sachdev, Phys. Rev. B 42, 4568 (1990). 





















6 


[26] K. Harada, N. Kawashima, and M. Troyer, Phys. Rev. 
Lett. 90, 117203 (2003). 

[27] The differentiation erases a summation as 

5 D /3 and produces the de¬ 

nominator, N. The numerator, M, comes from the 
Leibniz’s product rule. 

[28] S. Todo, in Strongly Correlated Systems: Numerical 
Methods (Springer Series in Solid-State Sciences), edited 
by A. Avella and F. Mancini (Sprineer-Verlae, Berlin, 
2013) pp. 153-184. 

[29] B. B. Beard and U. J. Wiese, Phys. Rev. Lett. 77, 5130 
(1996). 

[30] S. Todo and K. Kato, Phys. Rev. Lett. 87, 047203 (2001). 

[31] H. T. Quan, Z. Song, X. F. Liu, P. Zanardi, and C. P. 
Sun, Phys. Rev. Lett. 96, 140604 (2006). 


[32] D. Schwandt, F. Alet, and S. Capponi, Phys. Rev. Lett. 
103, 170501 (2009). 

[33] M. Kolodrubetz, Phys. Rev. B 89, 045107 (2014). 

[34] S. Chandrasekharan and U.-J. Wiese, Phys. Rev. Lett. 
83, 3116 (1999). 

[35] B. Bauer, L. D. Carr, H. G. Evertz, A. Feiguin, 
J. Freire, S. Fuchs, L. Camper, J. Gukelberger, E. Gull, 
S. Guertler, A. Hehn, R. Igarashi, S. V. Isakov, 
D. Koop, P. N. Ma, P. Mates, H. Matsuo, O. Parcol- 
let, G. Pawlowski, J. D. Picon, L. Pollet, E. Santos, 
V. W. Scarola, U. Schollwock, C. Silva, B. Surer, S. Todo, 
S. Trebst, M. Troyer, M. L. Wall, P. Werner, and S. Wes- 
sel, J. Stat. Mech.: Theo. Exp. , P05001 (2011). 

[36] “http://alps.comp-phys.org/,”. 



